home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / DDPOLY.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  724b  |  29 lines
  1. PROCEDURE ddpoly(c: glcarray; nc: integer; x: real;
  2.        VAR pd: glpdarray; nd: integer);
  3. (* Programs using routine DDPOLY must define the types
  4. TYPE
  5.    glcarray = ARRAY [1..nc] OF integer;
  6.    glpdarray = ARRAY [1..nd] OF integer;
  7. in the main routine. *)
  8. VAR
  9.    nnd,j,i: integer;
  10.    cnst: real;
  11. BEGIN
  12.    pd[1] := c[nc];
  13.    FOR j := 2 TO nd DO BEGIN
  14.       pd[j] := 0.0
  15.    END;
  16.    FOR i := nc-1 DOWNTO 1 DO BEGIN
  17.       IF (nd < (nc+1-i)) THEN nnd := nd ELSE nnd := nc+1-i;
  18.       FOR j := nnd DOWNTO 2 DO BEGIN
  19.          pd[j] := pd[j]*x+pd[j-1]
  20.       END;
  21.       pd[1] := pd[1]*x+c[i]
  22.    END;
  23.    cnst := 2.0;
  24.    FOR i := 3 TO nd DO BEGIN
  25.       pd[i] := cnst*pd[i];
  26.       cnst := cnst*i
  27.    END
  28. END;
  29.